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Abstract: In massive gravity and in bimetric theories of gravity, two constraints are 
needed to eliminate the two phase-space degrees of freedom of the Boulware-Deser ghost. 
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ghost is completely absent from these non-linear massive gravity theories and from the 
corresponding bimetric theories. 
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1. Introduction and summary 

Generic theories of massive gravity and bimetric gravity contain a Boulware-Deser ghost 
mode [jl| which renders these theories unstable. To avoid this inconsistency, the equations 
of motion of a healthy, ghost-free massive gravity or bimetric theory must contain two 
constraints which eliminate both the ghost field and the momentum canonically conjugate 
to it. In |3[ d it was established that one such constraint, the Hamiltonian constraint, 
exists in a special class of massive gravity and bimetric theories. In [^, |3[ Q an associated 
secondary constraint was argued to arise as well. Here we rigorously prove the existence 
of the secondary constraint and obtain an explicit expression for it. 

The class of theories being considered are extensions of the massive gravity actions 
originally proposed in [||, || where they were formulated with a flat reference metric. In 
these works, the proposed actions were shown to be ghost-free in the decoupling limit 
and to possess a Hamiltonian constraint up to fourth order in perturbation theory. For 
complementary analyses, see J?], 0. For other recent, related work, see H 10, 11, 12, 13 
14, 15, |l^, [l7], [lgl| . For a recent review on massive gravity, sec [19]. 



The full non-linear analysis of the Hamiltonian constraint was performed in [Q, |3[ 
and is based on a reformulation of these models presented in |20| ], where they are also 
extended to an arbitrary reference metric. In Q a ghost-free bimetric theory based on 
these massive gravity theories was proposed in which the reference metric is promoted to a 
dynamical field. The existence of the Hamiltonian constraint for the bimetric theory was 
also established in jl. For recent work related to the bimetric theory, see [^, 22, p3fl . 
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More specifically, the analyses of [||, [3|, ||] showed that both the massive gravity and 
bimetric theories contain a constraint, say C = 0, that can be regarded as a generalization 
of the Hamiltonian constraint in General Relativity. A condition for consistency is the 
preservation of this constraint in time, dC/dt = 0. This condition will lead to a secondary 
constraint, provided the Poisson bracket of C with itself vanishes on the constraint surface, 
[C(x),C(y)} » 0. In [Q, |3], U this bracket was not computed, but was assumed to vanish 
considering that it vanishes in two very different limits of massive gravity: the zero mass 
limit, which is General Relativity, and the linearized limit, which is the linear massive 
Fierz-Pauli theory [24, 25 1. It was then argued that a non-trivial secondary constraint 
should exist in the non-linear theory, as an extension of the secondary constraint in the 
Fierz-Pauli theory. This led to the conclusion that the Boulware-Deser ghost is eliminated 
by the two constraints in the non-linear theory. 

However, it was recently argued in j2(| that in massive gravity, the Poisson bracket 
{C(x),C(y)| does not vanish at the non- linear level, hence the theory does not contain 
a secondary constraint. If true, this would imply that, even though the ghost field is 
eliminated by the Hamiltonian constraint, the momentum canonically conjugate to the 
ghost is not eliminated. The theory would thus have an odd dimensional phase space - an 
odd situation indeed. In this paper, we explicitly compute the Poisson bracket in massive 
gravity and show that {C(x),C(y)} « indeed holds. Then we explicitly compute the 
secondary constraint C(2). We also argue that no other associated constraints exist. This 
result is then generalized from massive gravity to bimetric gravity. This conclusively proves 
the complete absence of the Boulware-Deser ghost both in non-linear massive gravity, as 
well as in non-linear bimetric gravity. 

The apparent contradiction of this result with the work |^6| can be understood in the 
following way. The constraint analysis of |2(| is based on a straightforward counting of 
Lagrange multipliers and equations of motion. However, as the existence of the Hamiltonian 
constraint itself demonstrates, such a counting can be misleading. If the equations of motion 
are degenerate, then an equation which would otherwise determine a Lagrange multiplier 
instead becomes a constraint on the dynamical variables. The analysis of [26] is insensitive 
to this possibility, while here we show that this is indeed the case. 

The paper is organized as follows: In section 2 we review non-linear massive gravity 
and its Hamiltonian constraint. Section 3 proves the existence of the secondary constraint 
and obtains an explicit expression for it. The absence of higher constraints is argued here 
as well. Section 4 extends the results to bimetric gravity. 



2. Review of the Hamiltonian constraint in massive gravity 

In this section we review the ADM formulation of massive gravity and the derivation of 
the Hamiltonian constraint based on [§,§]. 

2.1 Massive gravity in ADM variables 

We start by reviewing the derivation of the Hamiltonian constraint of the massive gravity 
theories proposed in [|, as was carried out in [§, Here, we construct these theories 
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with respect to an arbitrary but non-dynamical reference metric fav- In the last section 
of this paper we will discuss the generalization of these arguments to a theory with a 
dynamical f^ u . 

The most general massive gravity actions can be written as [^(| 

(2.1) 



r , r 3 , 

S = M 2 V I cfix^tetj i? + 2m 2 ^/3 n e n (V5 rT 7) 

n=0 



where the square root of the matrix is defined such that \fg~ =T f \[g~ =T ] = 9^f\u- The 
e k{\/ 9~ l f) ar e elementary symmetric polynomials of the eigenvalues A n of the matrix 
\J <7 -1 /. For a 4 x 4 matrix X and using the notation [M] = tr M, they can be written as, 

eo(X) = 1, 
ei(X) = [X], 

e 2 (X) = I([X] 2 - [X 2 ]), (2.2) 

ea(X) = i([X] 3 -3[X][X 2 ] + 2[X 3 ]), 

e 4 (X) = ^([X] 4 - 6[X] 2 [X 2 ] + 3[X 2 ] 2 + 8[X][X 3 ] - 6[X 4 ]) , 

e fc (X) = for fc > 4. 

The /3 n correspond to four free parameters, two of which are the graviton mass and the 
cosmological constant. 

The physical content of these theories is most easily identified in the ADM formulation 
1 27] which is based on a 3 + 1 decomposition of the metric. Let ./V and N{ denote the lapse 
and shift functions of the metric g^ v while L and Lj denote the lapse and shift functions 
of the metric f^ u , so that 

9^u- N = (-<?°°r 1/2 , Ni = g i, 7«=0ii, ( 2 -3) 
fc: L = (-/ 00 )- 1 / 2 , L i = f 0i , % = fij. (2.4) 

In massive gravity fu V is non-dynamical. The lapse and shift variables N„ appear without 
time derivatives in Q2.1| ) and are thus non-dynamical as well. The components jij are 
dynamical and, along with the their canonically conjugate momenta n 13 , constitute 12 
phase-space degrees of freedom, i.e., six potentially propagating modes. Of these, five 
correspond to the massive graviton while the sixth one, if not eliminated by the constraint 
equations, is the Boulware-Deser ghost. 

In terms of these variables, the Lagrangian in (|2.1| ) becomes, 



3 

£ = 7r^d tlij + NR° + R,N l + 2m 2 7det^iV £ /3 n e n {^/g^f) , (2.5) 

n=0 



where, as in General Relativity. 



R° = + -r^= (VVj " -' , --/) , Ri = 2>/diW7y Vjfe(-^=) • (2.6) 
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The equations of motion derived from fl2.5|) contain a single constraint on the variables 
7ij and tt 13 . This is the Hamiltonian constraint which is necessary (though not sufficient) 
to remove the ghost. However, since the potential in fl2.5| ) is highly nonlinear in the lapse 
N, the existence of a Hamiltonian constraint is far from apparent. To make this constraint 
manifest, we have to work with the appropriate set of variables. We define these in the 
following subsection. 

2.2 The new shift-like variables 

The existence of the Hamiltonian constraint is due to the fact that the four equations 
of motion in fact depend on only three functions n l (for i = 1,2,3) of the four variables 
Nu |2[ [|]. Thus three combinations of these equations determine the n l and the remaining 
equation is the Hamiltonian constraint. This can be made explicit if we work with the new 
shiftdike variables n % which are related to the ADM variables of ( |2.3|) and ( |2.4|) through 
(see pi for details), 



V = (LFj + ND'j) 



ir 



(2.7) 



where N l = r y lJ 'Nj and L l = 3 f lJ Lj. The matrix D = D" 1 - is determined by the condition, 

y/x D = J^r 1 - Dnn T D T ) 3 f , {2.1 

where, 



•'' 1 »' l : ]IV) n ! . (2.9) 

In (|2.8| ) we have introduced the notation n for the column vector n l , and n T for its trans- 
pose. We will occasionally use this matrix notation when the meaning is clear. The matrix 
D has the following important property 

%kD k j = % k D\ . (2.10) 

This identity will be used often to simplify our calculations. Although D can be determined 
in terms of ri 1 , the explicit solution is not needed here. 

We can use ([2.7D to eliminate the shifts A^ in favor of n l . The resulting form of the 
Lagrangian ( p.5| ) is then linear in N. It can be compactly written as, 

C = 7r ij d aij -Ho + NC. (2.11) 

Here Hq and C stand for 

n 



-{Ln l + U)Ri- 2m 2 L v / det 7 U, 
C = R° + RiD)n J + 2m 2 v / det7F . 
The functions U and V appearing above are given by the following expressions, 
U = fayfc + P 2 \V^ 2 D\ + n { %D\n k 



(2.12) 
(2.13) 



fa ^(D\n i %D\n k -D* k n h %D\n l ) + ^ (D^Di-D^) , (2.14) 



V^Po+Piy/xD'i + hfo^ 



D\D 



?>D i i D j k D k j + 2D i p\D k i 



(2.15) 
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The massive gravity Lagrangian is now given entirely in terms of 7^ , tt^ , N and n\ It also 
involves the ADM parameters of ( |2.4| ) as non-dynamical fields. 

2.3 The Hamiltonian constraint 



In the Lagrangian ( |2.11| ), the four variables N and n l are non-propagating while the six 
modes described by (7ij,7r lJ ) potentially propagate and include the ghost field and its 
conjugate momentum. The Hamiltonian constraint is identified after imposing the n k 
equations of motion. These were derived in 0, H|. Here we outline the derivation. Varying 



C in (2.11) with respect to n gives the equations of motion, 



dC 

dn k 



m> +N dc_ 

dn k dn k 



0. 



(2.16) 



To determine the variation of the two terms in ( |2.16| ), it is helpful to use the following 
relations, derived from the expression (|2.8|), 



a 



dn k 




-ti(\/xD) 



1 



T 3/. 

n f 



dn k 




: ti(y/xDy 



-2n T3 fD 



d(Dn) 
dn k ' 

d(Dn) 
dn k 



(2.17) 



dn k 



: tv(y/xDY 



One then obtains the useful expressions, 

on 



dn k 



-LCi 



dC 

dn k 



C 



dn k 



(2.18) 



Note that both these variations are proportional to Cj, given by, 



d = Ri- 2m 2 v / det 



■ n 



7- 



fij 



(2.19) 



Thus the variation of C becomes, 



dC 

dn k 



Ci 



L51 + N 



dn k 



0. 



The matrix within the square brackets is the Jacobian of the transformation 
invertible. Hence the n l equations of motion are 

Ci(7,7r,n) = 0. 



(2.20) 
and is 

(2.21) 



These are independent of the lapse N and can in principle be solved to determine n l in 
terms of jij and 7r* J . So far, the explicit solution is known only for the minimal massive 



action corresponding to (3% = P$ = 0. But, as shown below, the equations ( 2.21 ) alone are 
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enough to demonstrate the absence of the ghost. Note that ( p. 21 ) implies the vanishing of 



the individual variations in (2.18) 



Now, varying the Lagrangian with respect to N gives, 

C( 7 ,vr,n) = 0. (2.22) 



On eliminating the n* using the solutions 1 of ( |2.21 ), this equation becomes a constraint 



on the 12 dynamical variables jij and n 13 , reducing the number of independent degrees of 
freedom to 11. This is the Hamiltonian constraint. The degree of freedom it eliminates is 
the ghost field. One more constraint is necessary to eliminate the momentum canonically 
conjugate to the ghost field. This is obtained in the next section. 

3. The secondary constraint 

In this section we show that the Hamiltonian constraint gives rise to a secondary constraint 
and we obtain its explicit expression. This new constraint eliminates the momentum canon- 
ically conjugate to the ghost. Thus there remain 10 degrees of freedom out of an original 
12. 

3.1 Existence of a secondary constraint 

In order to be consistent, the Hamiltonian constraint C = must be preserved under the 
time evolution of 7« and Tr tJ , that is, dC/dt = 0. This requirement can be implemented in 
the Hamiltonian formulation. From the Lagrangian ( |2.1l| ), the Hamiltonian is given by, 

H = J d 3 x{n -NC). (3.1) 

In this expression the ri 1 are determined by their equations of motion ( |2.21| ), while N is a 
Lagrange multiplier. The only dynamical variables are jij and ir 13 . 

In the Hamiltonian formulation, the consistency condition on the time evolution of the 
constraint becomes, 

j t C(x) = {C(x),H}=0, (3.2) 
where the Poisson bracket is defined as, 

{C(,),g}-/A(7 C M " ,-/ C( ? " J. (3.3) 
1 V 1 5 J \hmn(z) S7r mn {z) Sn mn (z) 5~f mn (z)J v 1 

The 5 denote total variations, including dependence on jij and ir lJ through the n k . Using 
(|3.1|), the consistency condition (|3.2D becomes, 



{C(x),H} = J d 3 y{C(x),H (y)}- J d 3 y N(y) {C(x),C(y)} = . (3.4) 

It is necessary for this condition to hold only on the constraint surface, i.e., when C = 0. 
When imposing the above consistency condition, two possibilities could arise: 



lr That this is always possible follows from the fact that when the n l equations of motion are satisfied we 
have dC/dn 1 = 0. Hence on the space of solutions, C depends only on 7^ and n lJ . 
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(C(x), C{y) } 96 0, where the symbol ~ is used for equalities that hold on the constraint 
surface, i.e., when C = 0. In this case ( [3.4| ) becomes an equation for N and does not 
act as a constraint on 7^ and vr 1 -?. This leads to a theory with an odd dimensional 
phase space, or equivalently, with 5.5 propagating modes. Such a possibility was 
raised in pq] and was argued to be the case in massive gravity. 



2. {C(x),C(y)} « 0. In this case, the condition ( fO| ) is independent of iV and thus acts 
as a secondary constraint on 7^ and 7T*- 7 , 

C (2) (aO = {C(z),i?o} = 0. (3.5) 

where, Hq = J d 3 y7-Lo(y). The secondary constraint eliminates the momentum canon- 
ically conjugate to the ghost field and leads to a theory with 5 propagating modes. 
The Boulware-Deser ghost mode is thus entirely absent. This is the scenario advo- 
cated in g 1 1. 

In this section we will explicitly show that massive gravity indeed corresponds to the second 
case given above and is, therefore, entirely free of the Boulware-Deser ghost. We will show 
in section four that this argument also holds for the theories of bimetric gravity based on 
these massive gravity actions 0. 

To prove the existence of the secondary constraint, we now show that, 

{C(x),C(y)}^0. (3.6) 

To evaluate the Poisson bracket, note that C depends on 7^ and if 3 explicitly as well as 
implicitly, through the n % . However, from ( p. 18 ) it follows that when the n l equations of 



motion are satisfied one has j4 = 0- Thus the most general variation with respect to the 
canonical variables is given by, 

5C 



sc sc 



5j n 



5ir mn , (3.7) 



where the derivatives are now evaluated at fixed n % . This means that the Poisson bracket 
( |3.6j ) can be evaluated by replacing the total variations 5 by the partial ones evaluated at 
fixed n % . From now on we only consider these types of variations. 



Now we can compute |C(x),C(y)} using the expression ( 2.13 ) for C. To simplify, note 
that -R and V depend, respectively, only on 7r and 7 and not on their derivatives. Thus 
the term {R°(x), y/det 7 V(y)} is proportional to the Dirac delta function 5 3 (x — y). The 
net contribution of this term and its conjugate, -{R°(y),^/detjV(x)}, to the bracket 
{C(x),C(y)} is therefore zero, due to the antisymmetry of the bracket. Taking this into 
account gives, 

{C(x),C(y)} = {R°(x),R°(y)} + {R, t (x), R^y)} D\n k \x) D\n l \y) 
+ {R°(x),R l (y)}D i k n k (y) - {R°(y), Ri(x)} D\n k {x) 

+ S mn {x) D\n\y) - S™{y) ^| D\n k {x) . (3.8) 
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Here, the quantity S mn stands for 



S mn (x) = Rj{x) 



"linn v"1mn 



(3.9) 



The expression (|3.8|) involves the standard Poisson brackets of General Relativity (see, 
for example, p8|), 



{R°(x),R°(y)} 



{R*( x ),R i (y)} = -R°(y)£;5 3 (x-y), 



R\ x )^5 3 {x-y)-R i {y)^5 3 (x-y) 

0f„.\ d x3, 



(3.10) 



{Riix^Rjiy)} 



R {x)^5 3 (x -y)- Ri(y)^j5 3 (x - y) 



These terms as well as the remaining terms in ([3^) involve derivatives of the delta function. 
In order to manipulate these unambiguously, we introduce localized smoothing functions 
f(x) and g(y), and define 



F 



d 3 xf(x)C(x), 



G 



d 3 yg(y)C(y) 



(3.11) 



Then, the Poisson bracket {F, G} is given by 

{F,G} = [d 3 x[ d 3 yf(x)g(y){C(x),C(y)} 



(3.12) 



We will first compute this bracket and then from it extract |C(x),C(y) }. 

To determine the last two terms in ( |3.8| ), we note that, on using the expression (| 
for Rj, for any vector field = ^ k Vk we have, 



3 5Rj(x) 
d d x- J 



5ir mn (y) 



v j (x) = - [V m v n (y) + V n v m (y)} . 



(3.13) 



Using ( p.lO| ) and ( 3.13j ) in expression (|3.8|), and then carrying out one of the integrals in 
( 3.12[ ), one gets, 



{F,G} 



d 3 x 



f K d i9 + fdi(g R°D i k n k ) + f R^^d^g D 



+ 2fS mn V m (g lnj D\n k ) - (f g) 



(3.14) 



where everything under the integral is now a function of x. Again, due to the antisymmetry 
of the right-hand-side under the interchange of / and g, the only terms that do not cancel 
are the ones where a derivative directly acts on / or g. These can be written as 



d 3 x{fd l g-gd l f)P\ 



where, 



{F,G} 



P i = (R° + R j D j k n k )D\n l + R i + 2S il % j D\n k 



(3.15) 



(3.16) 
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To extract the Poisson bracket \C{x),C{y) } from ( 3.15 ), we write g(x) = f d 3 y g{y) 5 3 (x—y) 
and perform an analogous rewriting of /. Then, comparing with ( gig) gives, 



{C(x),C(y)} 



P\x) -^-5\x -y)- P\y) ^-5 3 (x - y) 



(3.17) 



The expression for P % can now be simplified. This is the key step in our proof of the 
existence of the secondary constraint. First, in the expression for S mn ( |3.9| ), we eliminate 
Rj in favor of Cj using ( 2.1 9| ) . We then compute dV/d^ mn using the following identities 
that are derived from the relation (|2.8|), 







d 




tr(v^-D) 

tl{y/xD) 2 



-2[n T3 fD^-^f d ^ 



2 <97„ 



This gives, 



where, 



d(D i n^) 

S mn (x) = C t 3 + m 2 ^/d^i (V~f mn - V mn ) , 



+ ft ( %k - %j D\ + \ % k (D' 1 ) k ] (£>', D\ - D\ D h { 



Since Ci is zero by the n l equations of motion, this is simply 



(3.18) 



(3.19) 



7 jn . (3.20) 



S mn (x) = m 2 v / det^ (F 7 mn - V mn ) . 



Thus P % ( |3.16|) becomes, 



p i = (r° + RjD j k n k + 2m 2 v / det7V r ) - 2m 2 y/d^yV il 'n j D 



(3.21) 



(3.22) 



The expression within the parentheses is simply C ( 2.13 ). Also, from the expression for 
ymn ( ft. 20 ) and on using ( 2.10 ), it is easy to verify that the remaining terms give Ci (|2.19j ). 
Thus, ( |3.22 ) is in fact, P % = C D 1 ^ 1 +Ci^f h . Again, since C, is zero on the n % equations of 
motion, this is equivalent to 

P i = C& l n l . (3.23) 
Thus the Poisson bracket fl3,17| ) finally becomes 



{C(x),C(y)} 



C(x) D^n\x) —^(x -y)- C(y) £>>'(y) —^{x - y) 



dy { 



(3.24) 
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As C(x) = on the constraint surface, it is clear that, 

{C(i),C(ff)}wO. (3.25) 

This corresponds to case 2 discussed following equation (^) and proves the existence of a 
secondary constraint. 

3.2 Evaluation of the secondary constraint 

The consistency condition dC/dt = now acts as a secondary constraint given by (|3.5|), 



C (2) (x) = J dZy{C{x),Uv(y)} =0. (3.26) 



To compute the Poisson bracket, note again equation (2.18), 



™} = -LC>, £=ft^£l (3.27) 

Given that vanishes on the n % equations of motion, we can therefor consider variations of 
both C and T-Lq at fixed n % . Then, using the expressions for Uo and C in ( p^) and 
we obtain, 

{C(x),H (y)} = -{R°(x),R i (y)}(Ln i +L i )(y) - {D\n k ){x) {Riix)^^)}^ +V){y) 

where we have introduced the notation, 

ITe (3. 29 ) 
V det 7 07 mn 



^mn - g gj ven gj ) an( \ a g a i n reduces to the expression ( 3.21| ). 



The brackets in the first line of ( |3.28 ) are the brackets of General Relativity given in 



( 3.10 ). Under the integral of ( 3.26[) , the ordinary derivatives appearing in ( 3.1 0| ) can be 



consistently converted to covariant derivatives. In addition, using Q and (glial) , we find 

S§) = ( 7mn(x)7r " fc(x) - 27w(x) ) 63{x - y) ' (3 - 30) 

5Ri(x) 

-{Hm(x)V y n +-f in (x)V y m)5 3 (x -y) . (3.31) 



5ir mn (y) 



To evaluate the secondary constraint we perform the integration of ( |3.26| ) . The derivatives 
of 5 3 (x — y) that appear in the above expressions can be manipulated under the integral. 



- 10 - 



After performing various integrations by parts to move the derivatives from the delta func- 
tions and after dropping the corresponding boundary terms, we determine the secondary 
constraint 



C (2) « m 2 L (j mn 4 - 2ir mn ) U mn + 2m 2 L^/detj (V m U mn ) lm D\n k 
+ (RjD 1 ^ - 2m 2 y / det~fj jk V ki ) V,(W + LP) 



+ Vdet 



7 



(Ln* + L*). (3.32) 



Here we have dropped a term CVi(Ln l + L 4 ) as this vanishes on the constraint surface. 
That the entire secondary constraint does not somehow vanish on the constraint surface 
can be seen by considering the form of the first term in ( fOp which appears nowhere in 
either C or the equation of motion Cj = 0. 

Now the consistency condition C/ 2 ) = can be used to eliminate the component of it 13 
that is canonically conjugate to the ghost component of 7y, which was itself eliminated by 
the primary constraint C = 0. This leaves only five propagating modes, establishing the 
absence of the Boulware-Deser ghost in massive gravity. 

3.3 Absence of a tertiary constraint 

It is important to verify that no further constraints are generated by the secondary con- 
straint and thus no additional degrees of freedom are eliminated 2 . In other words, the 
consistency condition 

j t C {2) (x) = {C {2) (x),H}=0, (3.33) 

should determine N as a function of the remaining variables, rather than act as a constraint 
on the dynamical variables. For this to be the case, we must have both 

{C {2) {x),H (y)} 76 and {C (2) (x),C(y)} 96O. (3.34) 

For both conditions, it is sufficient to consider the nonlinear theory at lowest order in 7 and 
7r. The theories considered here reproduce the Fierz-Pauli Hamiltonian at lowest order, by 
construction. Both constraints C and C( 2 ) reduce to the Fierz-Pauli primary and secondary 
constraints at lowest order. Since the above two conditions are satisfied by the Fierz-Pauli 
theory, they will also be satisfied by the full nonlinear theory. Thus the above consistency 
condition ( |3.33 ) becomes an equation for the lapse N, rather than a tertiary constraint. 



4. Extension to bimetric gravity 

It is straightforward to extend the above analysis to the case of bimetric gravity in which 
the reference metric fn V becomes dynamical. The most general bimetric theory is given by 



2 There are instances in which we expect massive gravity to have fewer than five propagating modes. For 
instance, massive gravity in the background of de Sitter spacetime is known to have only four propagating 
modes when the Higuchi bound is saturated Such a scenario is describable by the framework 

presented here. However, we expect that in these cases, the additional mode is removed by an additional 
gauge symmetry and not by tertiary and quaternary constraints. 
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a, 



S = M 2 g j d A x^-detgR (3) + M) J d A x^- det / 

+2m 2 M 2 ff / ^V-detg^^ e n {^g^]) . (4.1) 

n=0 

We have introduced two different Planck masses for the / and g sectors of the theory and 
defined an effective mass, 

The summation in the above action now runs from n = to 4. The contribution of the 
new n = 4 term is 



V- det g fa 64 ( vVV) = /W- det 7 • ( 4 - 3 ) 

Let us rewrite this action in terms of ADM variables, using the same shift-like variables 
n l introduced in (|2.7|), in place of the shift N l . The Lagrangian becomes, 



C = M] 7r«$7tf + M| p«$ a /ij -H + NC. (4.4) 
Here Hq and C stand for 

H = -V (M 2 g R { : ] + M 2 r\ S) ) - L (M}R oU) + M^Rf + 2m 2 M e 2 ffV / det^^) , (4.5) 

C = M 2 g R° {9) + M 2 g R[ 9) D)n 3 + 2m 2 M 2 s det 7 F . (4.6) 



In order to account for the new term, {5±yJ— det /, we have introduced 



y/det jU' = y/det jU + fa y/det 3 / , (4.7) 

where J7 and V in the above expressions are defined as before in ( p. 14 ) and ( p. 15 ) respec- 
tively. 

The p 13 are the momentum canonically conjugate to the dynamical variables 3 fij. The 
lapse L and shift L l of f^ u remain non-dynamical in the bimetric theory. Thus the 7« , n t3 , 
3 fij and p 13 represent 24 potentially propagating phase-space degrees of freedom. We argue 
now that the primary constraint C given above generates a secondary constraint, as in 
the massive gravity case, eliminating two phase-space degrees of freedom. In addition, we 
will argue for the existence of eight additional constraints, related to the restored general 
coordinate invariance of the bimetric theory, that remove an additional eight phase-space 
degrees of freedom. As a result, the bimetric theory has 14 phase-space degrees of freedom, 
or seven propagating modes, consistent with the seven modes one expects for a massive 
spin-2 field and a massless spin-2 field. 

The bimetric constraint C in (|0]) is identical to that defined above for massive gravity 
( 2.13 ) (up to irrelevant overall factors), but now depends on the dynamical variables 3 fij, 
both explicitly and through the shift-like variables n l . However, C is independent of p 13 , 
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the momentum canonically to 3 fij. To see this, note that the n l equations of motion, 
determined from the Lagrangian ( |4,4| ), are 

I 



M 2 q R { : ] -2m 2 M e 2 ffV / det^ n/ ' J 



m 

2 W iV^ n ^ ml 1 ^ i m 



0. (4. 



These equations can be used to fix the n l in terms of the dynamical variables. In addition 
to being independent of N, L and L l , these equations are also independent of p %3 . Thus 
n l = ri l ( , yij, 3 fij,'K 13 ). As a result, the constraint C is independent of p lJ as well. 

So, when calculating the Poisson bracket {C(x),C(y)} for bimetric gravity, one can in 
fact consistently treat the metric as non-dynamical, since the bracket evaluated with 
respect to ( 3 fij,p lJ ) is zero. Thus the proof of existence of the secondary constraint in 
massive gravity holds for bimetric gravity as well. The primary constraint C commutes 
with itself on the constraint surface. Thus the consistency condition dC/dt = acts as a 
secondary constraint C(2) on the dynamical variables, 7^, ir tJ , 3 fij and p 13 . 



In addition to these two constraints, we note that the Lagrangian (4.4) is manifestly 
linear in the lapse L and shift U . Thus the variation of the action with respect to these 
variables produces four additional constraints on the dynamical variables. These four con- 
straints, along with the four general coordinate invariances of the bimetric theory and the 
primary and secondary constraints C and C(2)> reduce the initial 24 potentially dynami- 
cal phase-space variables to 14, consistent with the seven modes of a massive spin-2 field 
coupled to a massless spin-2 field. 
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